Low-Loss All-Optical Zeno Switch in a Microdisk Cavity Using EIT 
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We present theoretical results of a low-loss all-optical switch based on electromagnetically induced 
transparency and the classical Zeno effect in a microdisk resonator. We show that a control beam 
can modify the atomic absorption of the evanescent field which suppresses the cavity field buildup 
and alters the path of a weak signal beam. We predict more than 35 dB of switching contrast with 
less than 0.1 dB loss using just 2 fiW of control-beam power for signal beams with less than single 
photon intensities inside the cavity. 



I. INTRODUCTION 

Over the past few decades transistors and other com- 
puting components have dropped in size while simulta- 
neously increasing performance. However, power dissipa- 
tion is increasingly becoming a fundamental limitation to 
performance [H . All-optical switches and transistors seek 
to address this issue, while simultaneously pushing for- 
ward the technology needed to create all-optical quantum 
computing devices [2-6]. One of the fundamental issues 
limiting this technology is the strength of the nonlinear 
effects that couple the signal and control fields. Electro- 
magnetically induced transparency (EIT) 0, H, has been 
investigated as a resource for optical switches [9l-[l2| and 
quantum memories [T^-EBl due to its large nonlinearity, 
which is enhanced by coherent effects. We demonstrate 
how it can be used along with a microdisk resonator to 
create a high-speed low-loss all-optical switch. 

The quantum Zeno effect (QZE) [l7| is a process 
whereby frequent measurements of a quantum system 
can inhibit transitions. In the first experimental demon- 
stration, the QZE was shown to inhibit driven transi- 
tions between 9 Be + ground-state hyperfine levels [18j| , It 
has been demonstrated that a classical analogue to the 
QZE can be used to create an all-optical switch using 
two-photon absorption (TPA) in a resonant optical cav- 
ity [HI, [2(| ■ This switch consists of a four-port resonator 
evanescently coupled to Rubidium (Rb) vapor. The pres- 
ence of two input beams results in sufficiently strong TPA 
to suppress the cavity field build-up thereby altering the 
path of the beams due to inteference effects. All-optical 
Zeno switches have also been proposed using other non- 
linearities including Raman [21[ and inverse Raman in- 
duced loss [22T ]. 

In addition to classical switching applications, it was 
shown that the Zeno effect could be a used for quantum 
information purposes. With a sufficiently strong nonlin- 
earity, it could enable quantum logic gates [23, 24] with- 
out the need for a large number of ancilla photons. The 
ability of the QZE to inhibit transitions could be used 
instead of protective ancilla photons to ensure that error 



events are suppressed in photonic quantum gates [2514291 ] . 
Related work has studied the QZE as a tool for protecting 
entanglement [30( and its effect on entanglement [3l| in 
general. It has also been suggested that the QZE could 
be used to implement quantum logic gates using semi- 
conductor materials [321 ] . 

Here we report theoretical results and performance es- 
timates showing that EIT and Autler-Townes splitting 
[33j can also be used to implement an all-optical switch 
in a microdisk cavity based on the Zeno effect. We make 
use of single photon absorption (SPA) to suppress the res- 
onant field buildup in a cavity, and use a control beam 
to modulate the absorption by inducing EIT (we refer 
to the combined EIT and Autler-Townes splitting effect 
as simply EIT from here out). The benefit of this ap- 
proach is that on-resonant SPA has a higher absorption 
cross-section than known nonlinear processes, potentially 
enabling better switching results as long as SPA can be 
sufficiently reduced by EIT. In addition this approach 
allows the use of very weak signal beams, making it a 
good candidate to switch single-photon intensities with 
low loss. 

To understand how the switch operates, consider the 
operation of a four-port resonator, as shown in Fig. [TJ in 
the absence of any atomic interaction (see e.g. [34| and 
35] for detailed discussions of four-port resonators and 
nomenclature). On resonance, a weakly-coupled signal 
Ei n will be almost completely transferred to the drop- 
port as a result of destructive interference between the 
resonant build-up of the cavity field and the light in the 
input waveguide. Next, consider how this system would 
change with the addition of a strong loss mechanism in 
the cavity, such as shown in Fig. 1(a) The resonant 
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field buildup of the signal beam will be suppressed which 
dramatically reduces the destructive interference in the 
input waveguide and results in nearly all of the incident 
light exiting the through-port. Counter-intuitively, the 
presence of the strong loss mechanism will not dramati- 
cally increase the loss of the system but will instead alter 
the coupling condition of the cavity which changes the 
output path of the light. 

This is directly analogous to the suppression of proba- 
bility amplitudes via measurement in the quantum Zeno 
effect. This is made more clear when considering the case 
of single photons in the signal beam. In that case, a null 
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Figure 1. Schematic of the EIT based Zeno all-optical micro- 
resonator switch. 

(a) When only a single input beam is weakly coupled to the 
cavity, strong SPA inhibits field buildup causing the beam to 
bypass the resonator and exit via the through port Et- 

(b) When two beams are present, denoted by blue for the 
signal beam and red for the control beam, the control beam 
eliminates the evanescent coupling of the fields to the atoms 
surrounding the cavity through EIT. This reduces the loss 
present, allowing the signal beam to build in the resonator, 
exiting through the drop port Ed- 



absorption (measurement) event by the atoms surround- 
ing the cavity causes the photon to bypass the resonator. 
The stronger this potential absorption process, the more 
likely that a null absorption event occurs and the pho- 
ton bypasses the c avity . More information on this can be 
found in Refs. fHHj. 

Now we consider a loss mechanism based on a cascade 
EIT scheme using the 551/2, 5P 3 / 2 and 5D 5 / 2 states of 
Rb. We take the signal beam to be resonant with the 
5S]/2 — * 5P 3 / 2 transition near 780 nm and the control 



beam to be resonant with the 5P 3 /2 — > 



5D 5 /2 transi- 



tion near 776 nm. When the 776 nm EIT control beam 
is present, a transmission window is created in the 780 
nm single-photon absorption line. This allows the signal 
beam to build in the resonator and exit through the drop 



port, as shown schematically in Fig. 1(b) 



In this paper, we theoretically analyze the performance 
of such a device using high-fidelity numerical models and 
present detailed estimates of switching performance. Our 
results indicate that this scheme enables high-contrast, 
low-loss switching at timescales on the order of the total 
cavity relaxation time which is roughly ^100 picosec- 
onds for the devices under consideration in this paper. 



Figure 2. Simulated normalized electric field profile of the 
microdisk, showing the radial component of the mode used in 
the quantum calculation. The cutout is the along page (side) 
view of Fig. [T] The radius of the device is set such that the 
free spectral range is equal to 4 nm, allowing for simultaneous 
resonance at 780nm and 776nm, and the thickness of the disk 
is chosen such that 30% of the field is outside the resonator. 



II. THEORETICAL MODEL 



We model the transmission characteristics of the four- 
port microdisk shown in Fig. [1] evanescently coupled to 
Rubidium vapor. The device specifications have been 
chosen to be consistent with current fabrication capabil- 
ities. The design consists of an unclad, suspended Si 3 N4 
disk with a free spectral range equal to the 4 nm differ- 
ence between the 55i/2 — > 5P 3/ /2 Di line at 780 nm and 
the 5P 3 / 2 5P/ 5 / 2 line at 776 nm, allowing for simulta- 
neous resonance at 776 nm and 780 nm (for an example 
of a similar device operating at a higher wavelength see 
Ref. [20]). The thickness is chosen such that roughly 
30% of the electric field energy of the fundamental mode 
is outside the resonator, allowing considerable interaction 
with the Rubidium vapor surrounding the cavity. 

We estimate the field profile of the fundamental mode 
in the cavity using a fully vectorial 2-D axially symmetric 
weighted residual formulation of Maxwell's equations im- 
plemented in Comsol Multiphysics software (as shown in 
Fig. |2J|. The field profile is used to estimate the evanes- 
cent interaction of signal and control beams with an en- 
semble of three-level atoms by calculating the effective 
absorption coefficient. This calculated absorption coeffi- 
cient is combined with a classical model of a resonator 
to predict switching performance. The intent is to show 
that the presence of the control beam will induce EIT, 
modifying the absorption coefficient of the signal beam 
in the resonator. This causes a change in the overall cou- 
pling between the waveguide and the resonator resulting 
in switching. 
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A. Atomic Model 

The device under consideration is designed to support 
simultaneous cavity resonances at the 780 nm and 776 
nm spectral lines of Rubdium. To model this interaction 
between the cavity fields and the Rubidium vapor, we 
approximate the Rubidium atom as the four-level atomic 
system shown in Fig. [3l The fourth level, indicated on 
the right, is the decay channel from the excited state 
level 1 3) to the ground state level |1) through the 6P 3 / 2 
level. It is included in our atomic model only as a decay 
term. The Hamiltonian of this system is the standard 
three-level cascade model given by 



H=(-tgi\i)(2\-ti%m\ 

+ ftA s |l)(l|-7iA c |3)<3|. 



h.c. 



(1) 



We denote all variables associated with the 776 nm con- 
trol field and 780 nm signal field with the subscripts c 
and s respectively. The Rabi frequency is f2 = 2d- i£/H, 
where d is the dipole moment, £ is the slowly varying 
electric field amplitude, e is the polarization vector, and 
A is the detuning of the field from its associated transi- 
tion. This Hamiltonian, along with the addition of phc- 
nomenological decay terms, gives the following density 
matrix equations: 
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where T mn is the homogeneous decay term associated 
with the m — ¥ n transition, and "f m n is the trans- 
verse decay term associated with the off-diagonal ele- 
ments. We take 712 = r 12 /2, 713 = (r i3 + r 23 )/2, and 
723 = (ri2 + r 13 + T 2 3)/2. These off-diagonal decay 
rates assume that all decoherence is due to population 
decay, which is suitable for atomic vapors. The values 
were chosen to be physically consistent with a positive 
density matrix, which does not always hold for arbitrary 
decoherence terms [3(1 H?} • The branching ratio between 
the T23 and R| 3 decay channels was set to 0.65 : 0.35 as 
given in Ref. [381 ] . 

Many analytic solutions to Eqs. © and similar models 
have been analyzed l39l - l4Tj . However, to ensure accu- 
rate device performance estimates, we use numerical solu- 
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Figure 3. The three-level atom we use to model the resonant 
5 Si/ 2 5P3/2 and 5P3/2 —¥ 5D 5 / 2 transitions in Rubidium, 
including all decay channels. The dotted lines indicate the 
5P3/2 splitting induced by the upper beam. The state on the 
right is the 6P3/2 which provides a decay channel from the 
5-D 5 / 2 excited state to the 5Si/ 2 ground state. 



tions to the density matrix equations above. This allows 
us to solve Eqs. ©, with only a steady-state approxima- 
tion, which is consistent with out experiments using CW 
lasers. However we do not make any perturbation ap- 
proximations since we assume a high-Q resonator, which 
can develop intense fields. The steady-state approxima- 
tion can easily be lifted to model pulsed phenomena, how- 
ever we include it here in order to reduce the run-time of 
our numerical simulations. 

We include Doppler broadening in our simulations to 
account for thermal motion of the atoms. We do so by 
averaging the density matrix elements over the Doppler 
profile given by 



F(A) = 



1 



2tt<td 



(3) 



where A is the laser field detuning due to the thermal 
velocity of the atom, and o\d is the Doppler linewidth, 
given by the Maxwell-Boltzmann velocity distribution. 



B. Cavity — Waveguide Coupling Model 

We model the coupling of the waveguides to the cavity 
using the coupled-mode equations of Haus [34[ (see also 
[35l V This method is related to the E&M approach taken 
in 19|, and we get nearly identical results using either 
technique. We will label the ports as done in Fig. [T] 
The input port will be called "in" and corresponds to the 
upper left waveguide in the figures, and the two output 
ports will be labelled the "drop" port and "through" port 
for the lower left and upper right waveguides respectively. 

The equation governing the field amplitude inside the 
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resonator for a single input field is is given by 

da 1 . — 

— = ~iAa - -(kq + n e + m + K 2 )a + iy/Kis in , (4) 

where a is the field amplitude inside the resonator, A = 
(lu—uq) is the detuning of the field with frequency to from 
the cavity resonance of frequency luq, kq is the intrinsic 
cavity linewidth, n e is related to the absorption rate of 
the atoms, K\ and k 2 are the cavity - waveguide coupling 
rates for waveguides 1 and 2 respectively, and |si„| 2 is the 
input power. 

The steady state solution to Eq. (@| is 



2^/k~i 



2iA + (k + K e + Kl + K 2 ) 



(5) 



The power out the through port is related to the power in 
plus the component that couples out from the resonator 

. 2lA + (/CQ + K e - Kl + K 2 ) 

St = s in + ty/^ia = ■ -s in . 6) 

2lA + (K + K e + Ki + K 2 ) 

The signal power out the drop port is related to the power 
that couples out from the resonator, 



Sd = i-JK 2 a 



~2y/KlK 2 



2iA + (Kq + K e + Kl + K 2 ) 



(7) 



The through port and drop port transmission rates 
can now be easily calculated. They are simply T = 
|st| 2 /|s m | 2 and D = \s d \ 2 /\s m \ 2 which yield 



T = 
D = 



4A 2 + (k + K e - K X + K 2 ) 2 
4A 2 + (k + K e + Ki + K 2 ) 2 

AKlK 2 

4A 2 + (Kq +K e + Kl+ K 2 ) 2 ' 



(8a) 
(8b) 



The intrinsic quality factor is defined as Qq = ujq/ko 
and the external quality factor is Qext — wo/(^e + ki + 
k 2 ). By controlling the atomic absorption rate K e we 
can modify the external Q of the atom-resonator system 
causing the coupling conditions to change, resulting in 
the ability to switch between the drop port and through 
port. 



C. Cavity Field — Atomic Interaction Model 

The interaction of the evanescent fields to the Rubid- 
ium surrounding the cavity is the source of the added loss 
K e in the waveguide-resonator coupling model above. To 
give performance estimates for the switch, we calculate 
the average absorption coefficient of the signal beam in 
the resonator. 

To calculate k £ we first solve the atomic equations, in- 
cluding Doppler broadening, given in ^ for a particular 
value of the signal and control beam fields. The base ab- 
sorption coefficient from standard quantum-optics theory 



is given by: 



a(Q, s ,Q, c ) 



ANd 2 LU X 12(^ S ,^c 

heoc 



(9) 



where Xi2(^s, ^c) = Im[pi2(^ s , fi c )]/^s is related to the 
linear susceptibility and is dependent upon: signal and 
control beam Rabi frequencies, the density N of Rubid- 
ium atoms, the dipole moment d and the angular fre- 
quency u) of the 780 nm transition. Here we explicitly 
note the functional dependence of all terms that depend 
upon the Rabi frequencies (field amplitudes) for clarity. 

The small mode volume and high Q of the microdisk 
produce large intra-cavity intensities for even modest in- 
put power levels. This can induce large splitting of the 
intermediate level. However, the field distribution of the 
control beam outside the resonator falls off sharply as 
shown in Fig. [5J To capture this, we compute the density 
matrix elements for each Rabi frequency corresponding 
to the different field values, and then average over the 
field distributions of both the signal and control fields 
as shown in Eq. ([TUl) below. These solutions are then 
averaged and weighted by the normalized signal beam 
intensity, to give the average absorption coefficient for 
the signal beam. This average absorption coefficient a is 
given by 



a(n s ,n b )w{r)dV, 



(10) 



where the integral is taken to extend over the volume ex- 
ternal to the resonator, denoted by V e , and the weighting 
function 



w(r) 



J Vt \n s (r)\ 2 dV 



(11) 



scales with the signal beam intensity. 

The normalization integral in the denominator of the 
weighting function extends over the entire volume of inte- 
gration Vt, including the region interior to the resonator, 
as opposed to the integral in Eq. (fTU)) which is only over 
the region exterior to the resonator. This accounts for the 
fact that the field inside the resonator does not interact 
with the Rubidium thus reducing the overall strength of 
the interaction. This weighting also accounts for the fact 
that the scattering cross section is larger near regions of 
high signal beam intensity than regions of low intensity. 
We assume that both beams are contained in the same 
cavity mode for this calculation. 

We can relate the average absorption coefficient, a to 
the external cavity loss rate K e through the simple rela- 
tion K e = ca, where c is the speed of light in the cavity. 
Thus, by changing the atomic absorption coefficient a 
through EIT, we can modify the total cavity Q, resulting 
in changing coupling conditions and switching. 
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Table I. Parameters used in atomic simulation. 




Signiil Beam Detuning A (GHz) 



Figure 4. Average absorption coefficient, denoted by a in 
the text, of the signal beam with the control beam on (blue 
curve) and control beam off (red curve) , plotted as a function 
of the signal beam detuning, A s . With the control beam on, 
the standard Doppler broadened line is split. Because of the 
non-uniform control field intensity, the usual Autler-Townes 
splitting is replaced with the absorption profile shown in blue. 




(a)Through port transmission (b)Drop port transmission equal 
equal bandwidth bandwidth 




Signal Beam Detuning A iV.Hzl 



iignal Benin Detuning A iViHz) 



(c)Through port transmission (d)Drop port transmission equal 
equal contrast contrast 

Figure 5. Through and drop port transmission plots for the 
parameters given in Tables ITU and ITTT1 Each plot provides the 
transmission percentage as a function of the signal laser de- 
tuning from the cavity resonance. The red curve corresponds 
to the case where EIT is off, and the blue curve corresponds 
to the case when the EIT control beam is on. In figures (a) 
and (b) we use the parameters in Table [IT] to equalize the 
bandwidth in both ports. In figures (c) and (d) we use the 
parameters in Table lTlll where the coupling condition was cho- 
sen such that the contrast was equalized. 



III. RESULTS 

Using the methodology above, the average absorption 
coefficient of the signal beam is calculated and plotted in 
Fig. 01 with the control beam on and off (blue and red 



Variable 


Description 


Value 


Units 


Pc 


Control Beam Power 


2.0 


MW 


p, 


Signal Beam Power 


20 


fW 


A 


Effective Mode Area 


8.2 x 10~ 10 


cm 2 


N 


Rubidium Density 


5.0 x 10 12 


cm -3 


d s 


Lower dipole moment 


2.1 x 10~ 29 


m-C 


d c 


Upper dipole moment 


4.6 x 10~ 30 


m-C 


Tl2 


Lower decay rate 


6.0 


MHz 


T23 


Upper decay rate 


280 


kHz 


Tl3 


Upper decay rate 


150 


kHz 


<7D 


Doppler linewidth 


240 


MHz 



curves respectively). All relevant atomic parameters for 
this simulation are given in Table [TJ The spectral shape 
of the absorption coefficient is modified from the typical 
double-Lorentzian peaks associated with Autler-Townes 
splitting because of Doppler broadening and the trans- 
verse variation of control beam intensity in the evanes- 
cent field. This causes the signal beam to see a distri- 
bution of EIT windows rather than a single split line as 
would be the case in a uniform field. Very large Autler- 
Townes splitting on the order of GHz is possible, because 
the high-Q and small mode volume of the resonator can 
generate high-intensity fields with relatively modest in- 
put powers. 

Upon calculating the average absorption coefficient of 
the signal beam in the resonator, we estimate the trans- 
missio n of the through-port (Fig. 5(a) I and drop-port 
(Fig. 5(b) I with and without the control beam, using 
©• 



the formulas in Eq. (jH]). We take the waveguides to 
be strongly over-coupled to the cavity with K\ and k-i 
roughly two orders of magnitude larger than «o . The cou- 
pling rates were chosen to equalize the bandwidth on the 
through-port and drop-port. With this constraint, we es- 
timate the on-resonant switching contrast to be 50 dB in 
the through-port and 25 dB in the drop-port, along with 
only 0.5 dB and 0.02 dB of loss in each respectively. We 
define the bandwidth as the point at which the switching 
contrast reaches 20 dB, on Figs. |5(a) and 5(b) This 



gives a bandwidth of approximately 516 MHz for each 
output port. These results are tabulated in Table [TT1 

Switching performance trade-offs can be analyzed by 
varying the coupling rates between the waveguides and 
cavity. In the previous example, we chose to equalize the 
bandwidth in the through-port and drop-port. This re- 
quired strongly over-coupling the waveguides to the cav- 
ity. We can instead choose to equalize the contrast, and 
by extension the loss in each port. To do this, we lower 
the fixed waveguide-cavity coupling rates which could be 
done by increasing the distance between the waveguides 
and the cavity for example. These results are shown in 
Figs. 5(c) and |5(d)"j with the corresponding performance 
metrics given in Table IIII1 On cavity resonance we esti- 
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Table II. Switch performance results - equal bandwidth 



Variable 


Description 


Value 


Units 


On 


Intrinsic Quality Factor 


3.6 x 10 6 




K\ 


Waveguide 1 coupling 


26.7 


GHz 




Waveguide 2 couplng 


26.7 


GHz 




R psiilt^ 

_l_ l_, V_. o LLJL Mj kj 








Through port loss 


0.5 


dB 




Through port contrast 


50 


dB 




Drop port loss 


0.02 


dB 




Drop port contrast 


25 


dB 




Through port bandwidth 


500 


MHz 




Drop port bandwidth 


500 


MHz 



mate 38 dB switching contrast with only 0.1 dB loss in 
each port. The tradeoff for equalized contrast, and lower 
loss is a reduced bandwidth in the through-port of 330 
MHz, however the drop-port bandwidth is increased a 
comparable amount due to conservation of energy con- 
straints. 



resonator. Because of the large resonant build-up that 
the micro-resonators allow, even with modest input pow- 
ers, very high intensity control fields can be created which 
can create excellent switching contrast with low loss. In 
addition, this method is compatible with a very low input 
power signal. These properties suggest that such a device 
could be suitable for switching single-photon intensities, 
ideal for some quantum information applications. 



Table III. Switch performance results - equal contrast 



Variable 


Description 


Value 


Units 


Qo 


Intrinsic Quality Factor 


3.6 x 10 6 






Waveguide 1 coupling 


5.9 


GHz 




Waveguide 2 couplng 


5.9 


GHz 




Results 








Through port loss 


0.1 


dB 




Through port contrast 


38 


dB 




Drop port loss 


0.1 


dB 




Drop port contrast 


38 


dB 




Through port bandwidth 


170 


MHz 




Drop port bandwidth 


840 


MHz 



IV. CONCLUSIONS 

We have presented numerical simulations demonstrat- 
ing the effectiveness of using EIT together with the classi- 
cal Zcno effect to create an all-optical switch in a micro- 
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